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ABSTRACT 

Molecular hydrogen is now understood to be the main coolant of the primordial gas clouds leading 
to the formation of the very first stars and galaxies. The line emissions associated with molecular 
hydrogen should then be a good tracer of the matter distribution at the onset of reionization of 
the universe. Here we propose intensity mapping of H2 line emission in rest-frame mid-infrared 
wavelengths to map out the spatial distribution of gas at redshifts z > 10. We calculate the expected 
mean intensity and clustering power spectrum for several H2 lines. We find that the 0-0S(3) rotational 
line at a rest wavelength of 9.66 /im is the brightest line over the redshift range of 10 to 30 with an 
intensity of about 5 to 10 Jy/sr at z ~ 15. To reduce astrophysical and instrumental systematics, we 
propose the cross-correlation between multiple lines of the H2 rotational and vibrational line emission 
spectrum. Our estimates of the intensity can be used as a guidance in planning instruments for future 
mid-IR spectroscopy missions such as SPICA. 

Subject headings: cosmology: theory — diffuse radiation — intergalactic medium — large scale struc- 
ture of universe 



1. INTRODUCTION 

Existing cosmological observations show that the 
reionization history of the universe at > 6 is likely 
both complex and inhomogeneous (e.g. Haiman 2003; 
Choudhury & Ferrara 2006; Zaroubi 2012). While the 
polarization signal in the cosmic microwave background 
(CMB) anisotropy power spectrum constrains the to- 
tal optical depth to electron scattering and the exist- 
ing WMAP measurements suggest reionization happened 
around Zri = 11 (Komatsu et al. 2011), it is more likely 
that the reionization period was extended over a broad 
range of redshifts from 20 to 6. Moving beyond CMB, 
observations of the 21-cm spin- flip line of neutral hydro- 
gen are now pursued to study the spatial distribution of 
the matter content during the epoch of reionization (e.g., 
Madau et al. 1997; Loeb & Zaldarriaga 2004; Gnedin & 
Shaver 2004). Unlike CMB, 21 cm data are useful as they 
provide a tomographic view of the reionization (Furlan- 
etto et al. 2004; Santos et al. 2005). The anisotropy 
power spectrum of the 21 cm line emission is also a use- 
ful cosmological probe (Santos & Cooray 2006; McQuinn 
et al. 2006; Bowman et al. 2007; Mao et al. 2008). 

While the 21 cm signal is primarily tracing the neu- 
tral hydrogen content in the intergalactic medium dur- 
ing reionization, line emission associated with atomic and 
molecular lines are of interest to study the physical prop- 
erties within dark matter halos, such as gas cooling, star- 
formation, and the spatial distribution of first stars and 
galaxies. Motivated by various experimental possibilities 
we have studied the reionization signal associated with 
the CO (Gong et al. 2011), CH (Gong et al. 2012), and 
Lyman-a (Silva et al. 2012) lines. As the signal is sensi- 
tive to the metal abundance these atomic and molecular 
probes are more sensitive to the late stages of reioniza- 
tion, perhaps well into the epoch when the universe is 
close to full reionization and has a low 21 cm signal (Basu 
et al. 2004; Righi et al. 2008; Visbal & Loeb 2010; CariUi 
2011; Lidz et al. 2011). 



Although the end of reionization era can be effectively 
probed with HI, CO, CH, and Lyman-a, it would be 
also useful to have an additional probe of the onset of 
reionization at z > 10. Here we consider molecular hy- 
drogen and study the signal associated with rotational 
and vibrational lines in the mid-IR wavelengths. Molec- 
ular hydrogen has been invoked as a significant coolant 
of primordial gas leading to the formation of first stars 
and galaxies (e.g. Haiman 1999; Bromm & Larson 2004; 
Glover 2005; Glover 2012). While molecular hydrogen is 
easily destroyed in later stages of reionization, its pres- 
ence in the earliest epochs of the cosmological history 
can be probed with line emission experiments. 

This paper is organized as follows: in the next sec- 
tion, we outline the calculation related to the cooling rate 
of H2 rotational and vibrational lines. In Section 3 we 
present results on the H2 luminosity as a function of the 
halo mass, while in Section 4, we discuss the mean H2 in- 
tensity and clustering auto and cross power spectra. The 
cross power spectra between various lines are proposed as 
a way to eliminate the low-redshift contamination and in- 
crease the overall signal-to-noise ratio for detection. We 
discuss potential detectability in Section 5. We summa- 
rize our results and conclude in Section 6. We assume 
the fiat ACDM with = 0.27, fl^ = 0.046, cts = 0.81, 
TT-s — 0.96 and h = 0.71 for the calculation throughout 
the paper (Komatsu et al. 2011). 

2. H2 COOLING COEFFICIENTS 

The radiation emitted by H2 will be generated by 
the heating/cooling of the gas as the collapsing process 
evolves. Therefore, in order to calculate the H2 lumi- 
nosity, we first evaluate the cooling rate of the H2 ro- 
tational and vibrational lines for optical thin and op- 
tical thick media, respectively. For hydrogen density 
riH < lO^cm"^, the optical depth is thin for H2 emission 
hues. Following HoUenbach & McKee (1979), the cool- 
ing coefficient for the rotational and vibrational lines can 
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then be expressed as 
A^i:(H,H2) = 



Alte(H, 
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where Alxe(H,H2) is the cooHng coefficients of rota- 
tion or vibration at the local thermodynamic equilibrium 
(LTE) for the collisions with hydrogen atoms H or H2. 
Here nl'^ is the critical density of H or H2 to reach the 
LTE, and tih.Hs the local number density of H or H2. 

In the LTE we have Au\ = Cu\n where A^i is the Ein- 
stein coefficient, Cui is the coUisional de-excitation rate 
from upper to lower level, and n is the particle number 
density (HoUenbach & McKee 1979). Then the rotational 
and vibrational LTE cooling coefficient can be written as 
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where gj — 2J + 1 is the statistical weight, J denotes 
the total angular momentum quantum number of the ro- 
tational energy level, Aj is the Einstein coefficient for 
J ^ J' transition at the same vibrational energy level 
or between two vibrational energy levels, and AE is the 
energy difference between Ej and Eji . The wavenum- 
ber for each energy level and the calculation for Ej and 
wavelength are given in the Appendix. The values of Aj 
are taken from Turner et al. (1977). Note that we just 
consider two vibrational level transition v=0 and 1 in the 
following calculation. The number of H2 at the higher vi- 
brational levels are much smaller than that at v=0 and 1 
in our case, so the strength of these lines is much weaker 
than those from v=0 and 1. The ricr/n term in Eq.(IT]) 
can be approximated by (Hollenbach & McKee 1979) 
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where Aj^'^g is the low-density limit of the cooling coef- 
ficient, which can be obtained by replacing the Einstein 



coefficient in Eq.® by Cj'^^nn.Haj i-^ 
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Here Cj'^^ is the coUisional de-excitation coefficients 
with H or H2 for J ^ J' transition in the same vi- 
brational level or between v=j and v=j, which are es- 
timated by the fitting formula given in Hollenbach & 
McKee (1979) and Hollenbach & McKee (1989) (see Ap- 
pendix). Using Eq.lllj-Q we can estimate the cooling 
coefficient for a given H2 line. Note that in Hollenbach 
& McKee (1979) the fitting formulae of total cooling co- 
efficients for v=0 rotational and v=0, 1 and 2 vibrational 
lines are given. We denote those by A^*'"^^, A^^g^ and 

A^°!^J7 cases of the total local, LTE and low-density 

limit cooling coeffecients, respectively. Then ricr/n can 
be expressed by AL'^g^/A„°_^g^, where ricr denotes the 
critical density when all energy level transitions are in 
the LTE. 

For the hydrogen density nn > lO^cm"^, the optical 
depth is thick for the H2 emission, and we have to con- 
sider the absorption effect for the H2 cooling. Following 
Yoshida et al. (2006), we make use of the cooling effi- 
ciency, which is defined by fee = Athick/Athin, to evaluate 



the cooling coefficient Athick for the optical thick case. 
This reduction factor is derived from their simulations, 
and available for nn < 10^'*cm~^, which is well within 
the density ranges of our calculation. The fee is about 
0.02 when nn ~ lO^^cm""^ and increases to about 0.1 
when TT-H ~ lO^^cm"^. At densities below 10^°cm~'^ we 
have fee = 1 (Yoshida et al. 2006). 
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Fig. 1. — The optical-thin H2 cooling coefficient vs. temperature 
for H2-H and H2-H2 collisions at z = 15. The cooling coefficients 
of the rotational lines, O-OS(O) (red), 0-OS{l) (magenta), 0-0S(2) 
(green) and 0-0S(3) (blue), are shown in dashed curves, and the 
vibrational lines, l-OS(l) (red) and l-OQ(l) (blue), are in solid 
curves. 

In Fig. [H as an example, we show the cooling coef- 
ficients for H2-H and H2-H2 collisions as a function of 
temperature T. We assume the number density of hy- 
drogen atom H and molecular hydrogen H2 to be 10^ 
cm~^ here. The dashed curves show the cooling coeffi- 
cients of the rotational lines, which are in red (O-OS(O)), 
magenta (O-OS(l)), green (0-0S(2)) and blue (0-0S(3)). 
The solid curves are for two vibrational lines l-OS(l) in 
red and l-OQ(l) in blue, respectively. As can be seen, 
the rotational cooling dominates at the low temperature 
(T < 10"^ K) and the vibrational cooling dominates at 
high temperature (T > 10'^ K). Also, we find the AJj^ for 
H2-H2 collision is generally greater than that for H2-H 
colfision at T < 10^ K, while the AJj^ for H2-H2 coUision 
is less than that for H2-H collision in this temperature 
range. This indicates that, at low temperature, the to- 
tal A|j^ is mainly from H2-H2 collisions, and the total 
Ag^ is from H2-H collisions. At higher temperature with 
T > 10^ K, the cooling rates A'^^ and A^^ for both H2-H2 
and H2-H collisions are similar. 

3. H2 LUMINOSITY 

We now explore the H2 luminosity as a function of 
halo mass for rotational and vibrational lines. As dis- 
cussed in the last section, the H2 cooling coefficient Ahj 
is dependent on the local gas temperature and density 
of hydrogen and molecular hydrogen. To evaluate H2 
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Fig. 2. — Left: The gas temperature T and H2 fraction fn^ as functions of the gas density n, which are derived from the simulation 
results in Omukai (2001) and Yoshida et al. (2006). The uncertainties of the gas temperature and H2 fraction are shown in blue regions. 
Right: The density profile of the gas (blue solid line) and molecular hydrogen (blue dashed line) for the halo with M = 10^ Mq /h at 
2 = 15. The blue region shows the uncertainty of the H2 density profile estimated by the uncertainty of /hj in the left panel. 



luminosity vs. halo mass relation for molecular hydro- 
gen cooling within primordial dark matter halos, we first 
need to know the radial profile of the gas temperature 
and density within dark matter halos. 

Following the results from numerical simulations in- 
volving the formation of primordial molecular clouds 
(e.g. Omukai & Nishi 1998; Abel et al. 2000; Omukai 
2001; Yoshida et al. 2006; McGreer & Bryan 2008), we 
assume the gas density profile as 
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where we set rg = 1 pc and po is the normalization factor 
which is obtained by 



the results on T{n) and /h2(") from Omukai (2001) and 
Yoshida et al. (2006), which are available for n ~ 10^^ to 
lO^'^ cm"'^ as shown in the left panel of Fig. [2] The uncer- 
tainties of the gas temperature and H2 fraction are shown 
in blue regions. These uncertainties are evaluated based 
on the differences in the far-ultraviolet radiation back- 
ground from the first stars and quasars (Omukai 2001). 

As can be seen, the gas temperature does not 
monotonously increase with the gas density. For in- 
stance, it drops from T ^ 2000 to 200 K between n ~ 1 
and 10'* cm~'^ where molecular hydrogen density is ris- 
ing to ~ 10~^. This indicates that H2 cooling is 
starting to become important in this gas density range. 
At n ~ 10'* cm"'^, H2 cooling saturates and turns into 
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the cooling at the LTE. For n = 10 
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Here Mgas = {Vti,/^M)M is the gas mass in the virial 
radius of the halo with dark matter mass M, and the 
Tvir is the virial radius which is given by 
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Here Pvir(-z) — ^c{z)pcr{z) is the virial density, Pcr(^) = 
3i?^(z)/(87rG) is the critical density at z, H{z) is the 
Hubble parameter, and Ac(2;) = 187r^-|-82a; — 39a;^ where 
X = VIm{z) — 1. 

We then derive the number density of gas by n(r) = 
nuir) + niie{r). Here n}i{r) = /hp(?')/toh is the num- 
ber density of hydrogen, where /h = 0.739 is the hy- 
drogen mass fraction and mn is the mass of hydrogen 
atom. Similarly, nHo('") = (1 — /H)/9(f)/wHc is the num- 
ber density of helium, where rriHc is the mass of helium 
atom. Also, the temperature-density relation T{n) and 
the H2 fraction-density relation /h2("') = nn^/n can be 
derived from existing numerical simulations. Here we use 



most all of gas particles become molecular hydrogen due 
to the efficient H2 three-body reaction (Yoshida et al. 
2006), and we find /h^ ^ 0.5 by definition. At n ~ 10^° 
cm~^ with T ~ 10* K, H2 begins to dissociate and the 
fraction drops quickly to < 10^^ when n ~ lO^'^ 
cm-3 and T ~ 10^ K. 

Next, with the help of Eq. ([5]), we can evaluate the 
gas temperature and H2 fraction as a function of the 
halo radius, i.e. T(r) and /h2('')- Once these are estab- 
hshed, we can derive 71H2 {r), nii{r) and Ay^(r) which are 
needed for the H2 luminosity calculation. In the right 
panel of Fig. [21 we show the density profile of the gas 
and molecular hydrogen in blue solid and dashed lines 
for a dark matter halo with M = 10^ Mo/h at z = 15. 
The blue region shows the uncertainty of the H2 density 
profile which is derived by the uncertainty of in the 
left panel of Fig. [2l The gas density profile is a straight 
line with a slope of -2.2 as indicated by Eq. ([5]). On the 
other hand, the density profile of molecular hydrogen has 
a more complex shape which is dependent on the rela- 
tion between fn^ and gas density n. For the outer layer 
of the gas cloud (r > 10^^ pc), gas density is less than 
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10^ cm^^ and fn^ ^ 4 x 10""'. Here H2 density is much 
smaller than the gas density. For the inner region with 
2Q-3.5 < ^ < 10-2 pc, we find 10^ < n < 10^^ cm~^ 
and begins to rise up quickly with becoming 
close to the gas density. For the inner-most region at 
r < 10~^'^ pc, the gas density is greater than lO^'^ cm""^, 
and — 0.5 so that almost all of hydrogen end up 
forming molecular hydrogen. 

The luminosity of the H2 rotational or vibrational lines 
can then be estimated by 
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where n]i^{r) is the number density of the molecular hy- 
drogen that can emit at a given rotational or vibrational 
line at r. We first evaluate the total n-n^ at v=0 and 1 
states by condensing all the rotational levels at a given 
vibrational state to be a single vibrational level, rii — 
Tii-i exp[— Ai?i.i_i/(fcT)] where i = 1. Here go = gi — ^ 
for v=0 and 1, respectively, and AEio/k ~ 5860 K (Hol- 
lenbach & McKee 1979). Then we estimate nn^ for a 
given rotational energy level J in a vibrational level i by 
nj = nji [gj/gj') exp[—A£' j^j//(fcT)]. The fractions of 
the ortho and para states of total tt-Hs are assumed to be 
0.75 and 0.25, respectively, in our calculation. 
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Fig. 3. — H2 line luminosity vs. dark matter halo mass M at 
z = 15. To avoid crowding we select the first eight strongest lines 
to show here. In our calculations we find that the 0-0S(3) is the 
most luminous line, while lines such as 0-0S{5), 1-0S{1), l-OQ(l) 
and 1-00(3), are comparable for high halo masses. 

In Fig. m we show H2 luminosity of several lines as 
a function of the halo mass at z = 15. In these lines, 
we find that the rotational line 0-0S(3) at a rest-frame 
wavelength of 9.66 /.tm is the most luminous one. Other 
lines, such as 0-0S(5), l-OS(l), l-OQ(l) and 1-00(3), are 
also strong for halos with high mass (see also Table[T|) . As 
can be seen, for low halo masses with M < 10^ /h, the 
rotational lines are stronger than the vibrational lines. 



This is caused by the fact that the halos with low masses 
have lower mean gas temperature than the massive halos, 
and the rotational cooling is stronger than the vibrational 
cooling in such halos, as indicated by Fig[TJ 

4. H2 INTENSITY AND POWER SPECTRUM 

Given the relation between H2 luminosity and the dark 
matter halo mass, the mean intensity of the H2 lines can 
be expressed as (Visbal & Loeb 2010; Gong et al. 2011) 
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where we choose Mmin ^ 
mass function (Sheth & 
Xu^il + zy/H{z) when 
is the wavelength of 
results are not strongly 
minimum halo mass. If 
mass to the level of 10^ 
present here decrease by 



= 10 Mq/H, dn/dM is the halo 
Tormen 1999), y{z) = dx/dv = 
X is the comoving distance and 
H2 lines in the rest frame. Our 
sensitive to the exact value of 
we increase the minimum halo 
Mq//i, the mean intensity we 
a factor of ~ 2 for all H2 lines. 




Fig. 4. — The mean intensity of H2 lines as a function of redshift 
z. The blue region is the uncertainty of the intensity of 0-08(3) line 
which is estimated from uncertainties on the gas temperature and 
shown on the left panel of Fig[2] We find the 0-0S(3) line is the 
most luminous line for 10 < z < 30, and the slopes of the relations 
for the vibrational lines are generally steeper than the rotational 
lines. Note that for 5 < z < 10 we do not consider the dissociation 
effect of the molecular hydrogen by Pop II and Pop III stars. 

In Fig. m we show the mean intensity of the eight 
strongest H2 lines as a function of redshift z. The un- 
certainty in the intensity of the 0-0S(3) line is shown 
with the shaded blue region, which is derived from un- 
certainties in the gas temperature and in the left 
panel of Figl^l We find that the mean intensity of the 
0-0S(3) rotational line is the strongest for 10 < z < 30. 
This is because the 0-0S(3) line is the most luminous line 
for low-mass halos, which have a higher number density 
and dominate the halo distribution at z ~ 15. Also, 
the slopes of the intensity-redshift relation for the vibra- 
tional lines are steeper than that of the rotational lines. 



5 




10' 



10° 



^_io- 



10-= 



10-2 10-1 IQO 

k [Mpc/h] 

Fig. 5. — Left: The H2 clustering auto power spectrum at 2 : 
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Eight brightest Hnes are selected to show here. The blue region is the 
uncertainty on the clustering power spectrum for the 0-0S(3) line and is derived from the uncertainties in the gas temperature and /h2- 
Right: The H2 clustering cross power spectrum at z = 15. Here we choose the strongest 0-0S{3) line to cross correlate with next five strong 
lines. The propose cross-correlation effectively eliminates the astrophysical line confusion from low-redshift sources and other mid-IR lines. 
The red region shows the uncertainty for 0-0S{3) x0-0S(5), which is also estimated from the uncertainties in the gas temperature and /h2- 



since they have steeper slopes for cooling coefficients with 
temperature as shown in Fig. [1] However, we find the 
difference in slopes to become smaller for the rotational 
lines as J is increased, indicating that the high-J rota- 
tional lines have similar slopes with cooling coefficient 
when compared to that of the vibrational lines. 

We note here that we do not consider the dissociation 
effect of the molecular hydrogen by Pop II and Pop III 
stars in our calculation. We expect the formation of these 
stars to be important at 2: < 10 and that there would be 
significant amount of H2 that should be be dissociated 
by the UV photons emitting from the first stars. Thus 
II2 emission could be suppressed significantly at z < 10. 
At z ~ 15, there should still be some dissociation but we 
ignore it to obtain a safe upper limit estimate on the ex- 
pected II2 intensity for experimental planning purposes. 

Next we can derive the clustering power spectrum of 
the H2 lines, writing the intensity as Ih2{^) — -^Hail + 
bn^Si^)]- Here &H2 is the average H2 clustering bias, 
which can be estimated from 



TABLE 1 

The wavelength, AJ = J — J' , spontaneous emission 

COEFFICIENT Aj, MEAN BIAS AND MEAN INTENSITY FOR THE H2 
ROTATIONAL AND VIBRATIONAL LINES AT 2: = 15. 
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where b{M, z) is the bias factor for dark matter halos 
with mass M at z (Sheth & Tormen 1999). The H2 
clustering auto power spectrum is then given by 



the spacial distribution of the primordial clouds, 
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where Ps5{k,z) is the matter power spectrum, which is 
obtained from a halo model (Cooray & Sheth 2002). At 
a high redshift as z = 15, the structure of matter distri- 
bution is extremely linear and the 2-halo term dominates 
the power spectrum. 

We can also estimate the shot-noise power spectrum 
for the II2 lines, which is caused by the discretization of 



In Table [1] we tabulate the rest-frame wavelength, 
A J = J — J', spontaneous emission coefficient Aj, mean 
bias and mean intensity for 13 rotational and 4 vibra- 
tional lines at z = 15. The uncertainties of the mean 
bias and intensity are evaluated by the uncertainty of 
the gas temperature and /jja from the simulations. We 
find that the mean intensity of the 0-0S(3) rotational line 
at a rest wavelength of 9.66 /xm is the strongest among 
these lines, with a value of around 6 Jy/sr and a range 
from 3 to 10 Jy/sr. The other rotational lines such as 0- 
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Fig. 6. — The auto power spectrum of the 0-0S(3) hne (left panel) and the cross power spectrum of the 0-0S(3) x0-0S(5) (right panel) 
at 2 = 12 with the errors estimated for a SPICA/BLISS-like and lOx better SPICA/BLISS-like surveys in each plot. The noise power 
spectrum and shot-noise power spectrum are shown in long-dashed and dotted lines, respectively. 



0S(5), 0-0S(4), O-OS(l) and 0-0S(2) are also bright with 
total intensities of ~ 4.3, 2.0, 1.5 and 1.3 Jy/sr, respec- 
tively, at z = 15. The vibrational lines l-OS(l), l-OQ(l) 
and 1-00(3) have low mean intensities at the level of 
0.83, 1.0 and 0.99 Jy/sr, respectively. The mean bias 
factor of these lines lies between 2.6 (for O-OS(O)) and 4.2 
(for O-OS(ll)), and the bias factors of the rotational lines 
at higher rotational energy level are higher than that at 
lower level. This is because the lines with high J are 
stronger at higher mass halos where the temperature is 
larger. 

In the left panel of Fig. [Sj the clustering auto power 
spectrum of eight H2 lines at z = 15 are shown. We find 
that the shot noise power spectrum Pghot is relatively 
small compared to the clustering power spectrum Pciusj 
and would not affect the Pcius at the scales of interest. 
This is easy to understand if we notice that the halo mass 
function is dominated by halos with low masses which are 
more abundant. 

We also calculate the cross correlation between two dif- 
ferent H2 lines. Such a cross-correlation will reduce the 
astrophysical contamination from the other sources, such 
as low-redshift emission lines from star-forming galaxies, 
including 63 /im [01] and 122 /xm [Nil], among others. 
At z ^ 15, the dominant rotational line 0-0S(3) would be 
observed at a wavelength of 155 /im. Such a line would be 
contaminated by, for example, z ~ 0.3 galaxies emitting 
[Nil]. Thus the auto power spectrum would be higher 
than what we have predicted given that the line intensi- 
ties of [Nil] are higher than II2 lines. To avoid this as- 
trophysical line confusion we propose a cross-correlation 
between two rotational or rotational and vibrational lines 
of the H2 line emission spectrum. 

The cross clustering and shot-noise power spectrum for 
such two II2 lines i and j can be evaluated as 



and 



^ shot ~ 



dM 



dn L\ 



dM 4:TtDI inDl 



clus 



(13) 



(14) 

respectively. From these equations, we find that the the 
cross power spectrum should have a similar magnitude 
to the auto power spectrum. The clustering cross power 
spectra Pfj^^ for several H2 lines at z = 15 are shown 
in the right panel of Fig. [5l We choose the strongest 
0-0S(3) line to cross correlate with the other 5 bright 
lines, i.e. O-OS(l), 0-0S(2), 0-0S(4), 0-0S(5) and l-OQ(l). 
We find that the 0-0S(3)x 0-0S(5) is the largest cross 
power spectrum since they are brightest two lines. At 
z 15, then we would be cross-correlating the wave- 
length regimes around 110 and 155 /im. A search for 
mid-IR lines revealed no astrophysical confusions from 
low redshifts that overlap in these two wavelengths at 
the same redshift. Thus, while low redshift lines will 
easily dominate the auto power spectra of H2 lines, the 
cross power spectrum will be independent of the low- 
redshift confusions. In addition to reducing the astro- 
physical confusions, the cross power spectra also have the 
advantage that it can minimize instrumental systematics 
and noise, depending on the exact design of an experi- 
ment. 

5. DETECTABILITY 

In this section we investigate the possibility to detect 
these lines based on current or future instruments. We 
assume a SPICA-likeQ survey with 3.5 m aperture diam- 
eter, 0.1 deg^ survey area, 10 GHz band width, i?=:700 
frequency resolution, 100 spectrometers, and 250 hours 
total integration time and noise per detector o-pix — 
10^ Jy-\/s/sr at 100 /xm. Such an instrument corresponds 
to the latest design of the mid-IR spectrometer, BLISS, 
from SPICA (Bradford et al. 2010). 

^ http:/ /sci.esa.int/science-e/www/area/index.cfm?fareaid=105 
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In Fig. ini we show the errors of the auto power spec- 
trum of 0-0S(3) hne and cross power spectrum of 0- 
0S(3)x0-0S(5) at z=12 for two cases, a SPICA/BLISS- 
Hke and an experiment with 10 x better sensitivity than 
with the current design of SPIC A/BLISS with CTpix = 
10^ Jy-\/s/sr. The noise power spectrum from the in- 
strument and shot-noise power spectrum caused by the 
discrete distribution of the gas clouds are also shown in 
long-dashed and dotted lines, respectively. We estimate 
the noise power spectrum and the errors by the same 
method described in Gong et al. (2012). We find the 
signal to noise ratio S/N is 0.2 and 5.2 for the auto power 
spectrum of 0-0S(3) line in the two cases, and S/N = 0.1 
and 4.5 for the cross power spectrum of 0-0S(3) x0-0S(5). 
This indicates that the current version of SPIC A/BLISS 
does not have the sensitivity to measure the intensity 
fluctuation of the H2 lines at a redshift around z = 12. 
We find that the noise requirements suggest an instru- 
ment that is roughly 10 times better in detector noise 
than current SPICA/BLISS for a reliable detection. 



T 1 1 1 1 1 1 1 1 Tj 1 1- r 

lOx SPICA (lff-250 hr) | ■ 

SPICA-like (1(7-250 hr) \ j 

I ! 




-28 -26 -24 -22 



log,„(flux [W/m=]) 

Fig. 7. — The number counts of H2 sources for 0-0S{3) line, 
per deg^ with flux greater than a given value for different red- 
shift ranges. The uncertainty for 2 = 10 to 15 is shown in blue 
region, which is estimated by the uncertainties of the gas density 
and /h2- We also show the flux limit of a SPICA/BLISS pencil- 
beam survey and another with 10 X better instrumental sensitivity, 
for a la detection with 250 hours of integration time. 

In addition to measuring the intensity fluctuations we 
also explore the detection of the II2 point sources at high 
redshifts. In Fig. [71 we estimate the number of the H2 
sources, for 0-0S(3) line per deg^ with flux greater than 
a given value for three redshift ranges 10 < z < 15, 
15 < z < 20 and 20 < z < 25. The uncertainty for 
10<z<15is shown as an example which is derived 
from the uncertainties of the simulations. The flux lim- 
its of a pencil-beam survey with a SPICA/BLISS-like 
instrument and a lOx better SPICA/BLISS surveys for 
la detection with 250 hours of total integration time are 
also shown in vertical dash-dotted lines. We find it is 
hard to detect the II2 sources even for the redshift range 



10 < z < 15 using the SPICA/BLISS-like experiment. 
The number counts of the H2 sources atl0<z<15is 
around 10"^ per deg^ for the SPICA/BLISS-like survey. 
For the instrument that is 10 times better in sensitivity 
than SPICA, we find that, in this first estimate of II2 
counts, that we can aim to get about 10 sources per deg^ 
at 10 < z < 15. 

6. DISCUSSION AND CONCLUSION 

In this paper, we propose intensity mapping of H2 rota- 
tional and vibrational lines to detect the primordial gas 
distribution at large scales during the pre-reionization 
epochs at z > 10. At such high redshifts the molecu- 
lar hydrogen takes the role of main coolant that leads to 
the formation of first stars and galaxies and the detec- 
tion of II2 power spectrum can reveal details about the 
halo mass scales which first form stars and galaxies in 
the universe. 

We first estimate the cooling rates for both H2 rota- 
tional and vibrational lines with the help of fitting re- 
sults from Hollenbach & McKee (1979) and Hollenbach 
& McKee (1989). We find the rotational lines are domi- 
nant at low gas temperature, while the vibrational lines 
are stronger at high gas temperature. Also, the slope of 
the cooling coefiicient-temperature relation for the vibra- 
tional lines is steeper than that for the rotational lines. 
We then derive the gas number density, temperature and 
II2 fraction as functions of the halo radius and estimate 
the relation of the II2 luminosity and halo mass. 

Next we calculate the mean intensity for several H2 
lines at different redshifts and find the 0-0S(3) is the 
brightest line for 5 < z < 30 (~ 5.9 Jy/sr at z=15). 
Note that we do not consider the dissociation effects of 
the II2 by the Pop III and Pop II stars at 5 < z < 10 that 
could suppress the H2 emission significantly in this red- 
shift range. Finally, we evaluate the clustering and shot- 
noise of auto and cross power spectrum for the H2 lines 
at z=15. We find the 0-0S(3) xO-OS(5) is the strongest 
cross power spectrum at z = 15. We propose such a 
cross-power spectrum for an experimental measurement 
as it has the advantage that it can minimize astrophysical 
line confusion from low-redshuft galaxies. 

In order to consider potential detection of these mid-IR 
molecular lines we evaluate the errors of the H2 auto and 
cross power spectrum at z = 12 for a SPICA/BLISS-like 
and a design that is 10 x better than the current instru- 
mental parameters. We find the S/N for the z = 12 cross 
power spectrum detection is around 0.1 and 5 for these 
two experiments. We also estimate the detectability of 
II2 point sources over different redshift ranges. We find 
a SPICA/BLISS-like instrument is not able to detect II2 
sources for z > 10, but an instrument with 10 x better 
sensitivity than SPICA/BLISS should be able to detect 
about 10 sources per deg^ for 10 < z < 15. We encourage 
further work on this topic to fully account for dissocia- 
tion as stars and galaxies form and additional formation 
mechanisms, such as shock heating, that results in H2 
line emission from low-redshift galaxies. 

This work was supported by NSF CAREER AST- 
0645427. MGS acknowledges support from FCT- 
Portugal under grant PTDC/FIS/100170/2008. We 
thank Matt Bradford for helpful discussions and ques- 
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TABLE 2 

The wavenumbers in cm^-'- of H2 energy levels for J=0,...,13 at v=0 and 1. 



7 

J 


v=0 


v=l 





0.00 


4161.14 


1 


118.50 


4273.75 


2 


354.35 


4497.82 


3 


705.54 


4831.41 


4 


1168.78 


5271.36 


5 


1740.21 


5813.95 


6 


2414.76 


6454.28 


7 


3187.57 


7187.44 


8 


4051.73 


8007.77 


9 


5001.97 


8908.28 


10 


6030.81 


9883.79 


11 


7132.03 


10927.12 


12 


8298.61 


12031.44 


13 


9523.82 


13191.06 



tions that motivated this paper. 
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APPENDIX 



Here we Ust the fitting formulae of the collisional de-excitation coefficients Cj '^^ for both H2-H and H2-H2 collision 
from HoUenbach & McKee (1979) and HoUenbach & McKee (1989). For the rotational cooling in v=0, we have 



0.33 + 0.9 exp 



C^^T) = (3.3 X 10-12 -t- 6.6 X IQ-i^Tg) J 0.276 exp 



J -3.5 
0.9 

J 
3l8 



1.7 



Q —1 
cm s , 



Q _ 1 
cm s , 



where T3 = T/1000 K, and T is the gas temperature. For the vibrational cooling between v=l and 0, we have 

Cfo(T) = 1.0 X 10^12 exph(1000/r)] cm^s^S 
Cfo'(r) = 1.4 X 10^12 r"-5exp-[18100/(T+ 1200)] cm^s"!. 



(1) 
(2) 

(3) 
(4) 



The wavenumbers fc = 1/A of H2 energy levels for J=0,...,13 at v=0 and 1 from Dabrowski (1984) are listed in Table 
[21 We can get the energy for each level by Ej = hpck, where hp is the Planck constant and c is the speed of light. 
The wavelengths of the H2 line hence could be derived by Ahj — hpc/{Ej — Eji). 



